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ABSTRACT 



We present examples of fits of fermion mass data using the 10 — 120 — 126 FM 
Higgs system in Susy SO(10) GUTs that follow the scenario P |2j in which the 120 
-plet collaborates with the 10-plet to fit the charged fermion masses while small 126- 
plet couplings enhance the Type I seesaw neutrino masses to viable values and make 
the fit to light fermions accurate. Restricting ourselves to the CP conserving case we 
use a linear perturbative technique to obtain accurate charged fermion fits (xfot < -2) 
to all the charged fermion masses and angles. The resulting fits imply Type I neutrino 
masses that are generically 10 2 — 10 3 times larger than those obtained in the 10 — 126 
scenario precisely because of the small 126 coupling. Thus the difficulty of obtaining 
sufficiently large neutrino masses in the context of these Next to Minimal Susy GUTs 
is essentially removed. The remaining free parameters allow one to obtain the correct 
ratio of neutrino mass squared splitting and a large(6^ MAr5 ) and sma mixing 
angle. The 9f 2 M is however small and this indicates that -as in the Type 1 10 — 126 
case - a fully realistic fit to the lepton mixing data also requires CP violation. 



1 Introduction 

The discovery of neutrino oscillations in 1998, followed by the measurement of in- 
creasingly accurate values of the neutrino mass squared differences and mixing an- 
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gles has driven an intense wave of research into supersymmetric seesaw mechanisms, 
Left-right symmetric models and particularly all aspects of the structure of renor- 
malizable Susy SO(10) GUTs. In particular the minimal Supersymmetric Grand 
Unified Theory (MSGUT) jSHHEj namely a Supersymmetric SO(10) GUT based on 
a 210 — 10 — 126 — 126 Higgs system has received a great deal of attention over 
the last seven years. Detailed calculations of every aspect of its spectrajHl 13 El E|, 
couplings |B1 IE] and their evolution [HI E] have been published by several groups. Be- 
sides the question of gauge unification the central problem of the second phase of the 
Grand Unification program namely that of accounting for the observed fermion mass 
spectrum and mixing structure, for both charged fermions and neutrinos, using the 
mass relations dictated by the simplest possibly realistic 10 — 126 Higgs structure 
[TU] . has inspired a great deal of detailed analysis of the fitting problem [T2^ ITS! IT^ Ho] . 

This so called Babu-Mohapatra(BM) program had, until last year, met with uni- 
form success in accounting for the observed charged fermion spectrum, neutrino mass 
splitting and mixing using the generic structure of the fermion mass matrices dic- 
tated by the choice of Higgs (10 — 126) and the Type lb and Type II |20] Seesaw 
mechanisms naturally present in SO(10) models. Successful fits of all the fermion 
mass data available including the neutrino mixing data [22] using Type I , Type II 
and mixed seesaw mechanisms were obtained[TTl[T21[Tni EJUS], ft was > however, just 
assumed that the overall scale and the relative magnitudes of the Type I and Type 
II seesaw masses required for the fits would be realizable in appropriate Susy GUTs. 
The only such available theory, where in fact the calculations of spectra (HI EH III Ej 
and couplings [HJ El HE] necessary to perform precisely such a test had almost con- 
temporaneously been completed, was the MSGUT itself. The very first survey [TH] 
revealed that there were serious difficulties in obtaining Type II over Type I domi- 
nance as well as in obtaining large enough Type I neutrino masses. Using a convenient 
parametrization of the MSGUT spectra and couplings in terms of the single "fast" 
parameter (x) which controlled MSGUT symmetry breaking [3 El El HZ] a complete 
proof was then given 1 of the failure of the Seesaw mechanism in the context of the 
MSGUT. The nature of the obstruction uncovered by us also led us to suggest a very 
natural scenario that might deal with the problem in the MSGUT extended by adding 
the 120 FM Higgs representation(Next to Minimal Susy GUT or NMSGUT). The 
120 had been optimistically excluded in the Babu-Mohapatra program to foster the 
predictivity of the MSGUT in the neutrino sector. In this alternative scenariopQ the 
120 -plet collaborates with the 10 -plet to fit the dominant charged fermion masses. 
The small 126 -plet couplings give appreciable contributions only to light charged 
fermion masses and enhance the Type I seesaw masses to viable values since the 
Type one seesaw masses are inversely proportional to these masses. To begin anal- 
ysis of the fitting problem in this new and qualitatively different fitting scenario we 
analyzed [2] the 2-3 generation toy model of the dominant core of the com- 

plete hierarchical fermion mass system. We found that ansatz consistency requires 
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nib — m s = m r — m M at the GUT scale Mx and predicts near maximal (PMNS) 
mixing in the leptonic sector for central values of charged fermion parameters and 
for wide ranges of the other relevant parameters (righthanded neutrino masses and 
relative strength of contributions of the two doublet pairs from the 120-plet to the 
effective MSSM Higgs pair). 

In the current contribution we report on progress in extending the analysis of j2] 
to a 3 generation CP conserving but otherwise realistic NMSGUT model. We first 
briefly review the failure of the seesaw mechanism in the MSGUT and the results 
of our analysis of the 2-3 generation toy model. This is followed by a description 
of the fitting equations in the CP conserving 3 generation case and our procedure 
for obtaining accurate solutions of the 3 generation equations by an effectively linear 
perturbation expansion (in the Wolfenstein parameter e ~ sin#J 2 ~ y/sin #23 ~ -2) 
around the dominant 2-3 generation mass matrix structures. We then present ex- 
amples of solutions that are realistic in all respects except that of the omitted CP 
violation and the lack of two large mixing angles. We argue that since even in the 
10 — 126 Type I seesaw a solution with two large angles and the observed small ratio 
of mass splittings |22j could be found only when CP violation was also included |lH lTo] 
our results should be construed positively. The study of the CP violating case is now 
in progress [22] ■ Thus the main obstruction of too small Type I masses is already over- 
come even in the CP conserving model but the fit for the CP violating case remains 
to be done. 

2 Seesaw Failure in the MSGUT 

The Type I and Type II seesaw Majorana masses of the light neutrinos in the MSGUT 
are [I] : 

(1.70 x 10~ 3 eV) h sin/3 
(1.70 x 10~ 3 eU) F H f sin/3 

(h-3f)r\h-3f) 

where v = 174 GeV and h, f are the Yukawa coupling matrices of 10, 126 to the 
16 plets containing fermion families and j3 is the MSSM Higgs doublet mixing angle. 
The functions Fj, Fjj are completely specified by in the MSGUT and explicit forms 
may be found in El UJ Ej . 

In typical BM-Type II fits |13( ITU ITo] the maximal value of / eigenvalues is ~ 10~ 2 
while the corresponding values for h are about 10 2 times larger. As a result n max ~ 
10 2 . This implies that R = Fj/Fji < 10~ 3 in order that the pure BM-Type II not be 
overwhelmed by the BM-Type I values it implies. Such R values were shown jTHJ 
to be un-achievable anywhere in parameter space of the MSGUT (while preserving 
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baryon stability, perturbativity etc). Furthermore in the BM-Type I fit of [TT) ITo] 
the requirement of large neutrino mixing yields, typically, n max ~ 5f max ~ .5 for the 
maximal eigenvalues of n. Thus values of Fj ~ 100 are required in order to reach 
realistic values of seesaw masses for the heaviest neutrino. We demonstrated [TBI UJ 
that such values are not achievable anywhere over the complex x plane except where 
they also violate some aspect of successful unification. Recently another group has 
verified our results [T5] 



3 The new 10 - 120 - 126 scenario 

To resolve the difficulty with the overall neutrino mass sscale we proposed^! |2] that 
the 126 couplings be reduced far below the level where they are important for 2-3 
generation masses while introducing a 120 plet to do the work of charged fermion 
mass fitting previously accomplished by 126 couplings almost comparable to those 
of the 10-plet. 

The Dirac masses in such GUTs are then generically given byjfH El E E] 
m u = v(h + f + g) 

rh v = v(h-3f + (r 5 -3)g) =v(h-3f + r' 5 g) 

m d = v(r l h + r 2 f + r 6 g) (1) 

rh l = v {rji — 3r 2 f + r 7 g) 

See |H1 El UJ Ej for the form of the coefficients r«. 

The right handed neutrino mass is M D = fa and the Type I seesaw formula is 

Ml = vr,h ; h = (h-3f-r' s g)f- 1 (h-Zf + r' s g) (2) 

where a = a is the GUT scale vev of the 126 while <%■, a,- refer to fractions 

of the MSSM doublets contributed by various doublets present in the GUT Higgs 
representations 13 El UH Uj • 

The essence of our proposal is that the 10, 120 multiplet Yukawas dominate the 
contributions of the 126 coupling in the charged fermion masses i.e (h,g) max » 
f. We shall work with the assummption that the fermion Yukawa couplings and 
coefficients are both real. The latter assumption is almost certainly not valid jU El El- 
However if a fit could be found with all complexity introduced via these coeffcients 
and the Yukawa couplings were kept real then the resulting NMSGUT would contain 
only 12 Yukawa coupling parameters i.e 3 less than in the MSGUT with complex 
Yukawas. Thus even with the new couplings due to the 120 the NMSGUT with only 
spontaneous CP violation would still be a strong contender as far as minimality of 
parameters is considered as decisive. 

The mass terms above must be matched to the renormalized mass matrices of 
the MSSM evaluated at the GUT scale. While doing so one must allow^H] for the 
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possibility that the fields of the GUT are only unitarily related to those of the MSSM 
at Mx- This introduces several unitary matrices into the fitting problem which, 
besides conventional ambiguities, specify just how the MSSM lies within the MSGUT 
once conventions are fixed. Although unphysical in the MSSM these matrices are of 
vital relevance for calculation of the exotic signatures (e.g B-L violation) of the GUT 
in question |16j. Thus when matching the charged fermion dirac mass matrices we get 

m u = V?D U Q 
m d = VjD d R 

m l = VfDiL (3) 

Where D Ut d,i are the charged fermion masses at Mx and V u ,Q,Vd,R = C^Q,L,Vi 
are arbitrary unitary matrices (C is the CKM matrix). These matrices must be 
fixed by convention -where allowed by conventional ambiguities - or determined by 
the fitting procedure in terms of the low energy data, or left as parameters to be 
determined by future experiments sensitive to degrees of freedom and couplings (e.g 
baryon violating couplings) that the low energy data is not. We shall make the SO(10) 
basis choice R = 1 in what follows. 

To put our equations in a form transparent enough to clearly separate the con- 
tributions of the 10, 120 (and -eventually- 126) plets[2"lj we write the matrices V u ^,i 
associated with the anti-fermion fields as new unitary matrices <& u ,d,i times the fermion 
field matrices Q,R,L. 

V d = ; V U = $ U Q ; V t = ^L (4) 

We then separate symmetric and antisymmetric parts for each charged fermion equa- 
tion 

Z = $ T D + D$ ; A = $ T D- D$ (5) 
Solving the mass formulae for the symmetric yukawa couplings one finds : 

h = ^-R T (3Z d + VZ l V T )R 

f = ^R T (Z d -VZ l V T )R (6) 
while the antisymmetric couplings are 

9 = 7T—R T AiR (7) 

2vr 7 

and where fj = l/(4rj). 

Then the remaining equations are the "sum rule" 

(3f x + f 2 ) Z d + (f ! - h)VZ{D T - C T Z U C = (8) 
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3.1 23 generation toy model 



We first briefly recapitulate the analysis of the 2 generation CP conserving toy model 
given in |2j. As explained, in this case the contribution of / to the charged fermion 
masses is taken to be negligible. Then using (say) just the d-type quark equations to 
solve for h, g one obtains for the rest : 



C*Z d & = nZ u ; Z d = VZ{D T ■ V = R*L T 

C*A d & = r 6 A u ■ r 7 A d = r§T>A{D T (9) 

In this 2 generation case the antisymmetric equations© serve only to fix the 
parameters tq,tj and thus play no further role. On the other hand since C,T> are 
orthogonal matrices it follows that r% = Tr(Z d ) /Tr(Z u ) so that we can write these 
two equations in the (dimensionless) form 



3 CZ d C T Z u _ 

1 TrZ d TrZ u 

a = z '- pz ' pT ^o no) 

TrZ d 1 ' 

In [2] we gave an analytic solution of the toy model charged fermion mass fitting 
relations. We parametrized the matrices $ ud i,C,T> as 



dj 



cos Xu,d,i sin Xu,d,i 
- sin Xu,d,i cos Xu,d,i 



CV = I COSXc ^ sm Xc,v \ 
V - sin Xc,v cos Xc,v J 

Then we found that the equations S± = yield (a = (^^Xtfe-"^) ) 



Xu Xd ^Xc 

csc 2xc 

tanxd = cot2x c (12) 

a 

Since a = 1 + (e 2 ), it is clear that the leading (~ e -2 ) contributions cancel leaving 
behind an 0(1) result. Similarly two of the three equations §2 = yield 



Xd 
Xi 



(Xi - Xd)/2 

.! T d COSXd 



Xi = ±Xi 



cos 



(13) 
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To leading order xi — ±Xd, this defines two branches of the solution which we call 
the (+) and (-) branches following the sign between xi an d Xi- The third 5*2 = 
equation yields the all important consistency condition 

d 3 - d 2 = Ad = Al = l 3 - l 2 (14) 

i.e 

m b (M x ) - m T (M x ) = m s (M x ) - m^M x ) (15) 

Using this solution we determined the 23 sector PMNS mixing angle as a function 
of the parameters r f 5 , p = f 3 / f 2 (fi are the eigenvalues of /) and then showed that 
the 23 sector PMNS mixing was near maximal over much of the parameter space of 
the theory. Thus the toy model provides a cartoon of the core of the 23 generation 
dominated fermion hierarchy that is satisfactory and coherent with our scenario in 
all respects. 

3.2 3 generation fitting formulae 

We now wish to generalize our 2 generation toy model to the 3 generation case, while 
keeping the CP conserving approximation for simplicity, to enquire how closely one 
may approach the actual fermion data. The resulting equations are however too com- 
plicated to solve analytically. Thus one must resort to some variety of numerical 
approximation technique or perturbation theory. Our approach is to consider the 
23 sector as the dominant "core" of the fermion mass hierarchy and thus to expand 
around this core by rewriting all the charged fermion mass parameters and mixing 
angles relative to the dominant normative magnitudes of the third generation by in- 
troducing a small parameter e ~ \fX c ~ sm ®\i ~ V^i3 ~ \J~d2Jd3 etc. The explicit 
analytic solution of the toy model given above coincides, order by order in pertur- 
bation theory , including the b — r = s — ji constraint , with the solution found by 
expanding all angles x m powers of e and solving the equations Si = = S 2 order 
by order in e. The b — r = s — /i unification constraint arises at order e 2 . In the 
three generation case the analytic solution is not available but the complete regular- 
ity observed in the perturbation theory in e and the extreme smallness of the first 
generation perturbations makes an expansion in e well motivated. Thus we shall cast 
the fitting equations in a form such that they reduce to the toy model equations up 
to order e 2 since /, d\, u±, l\ appear only at order e 3 or even higher. 

We continue with our CP conserving restriction to maintain tractability and 
parametrize each of the orthogonal matrices in the standard Kobayashi-Masakawa 
form for real unitary, i.e orthogonal, matrices : 

O = 23 ( X ) ■ OM) ■ O 12 (0) (16) 
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where Xi <$>■> ® are generic symbols for rotation angles in the 23, 13, 12 sectors. We 
shall always work in the SO(10) basis fixed by the condition R — 1. 
Now the antisymmetric equations are no longer empty but read 



\J ^i,d,l 

CA d C T - A u = 

A d =F VA{D T = (17) 
where 

A 2 = A 2 12 + A 2 13 + A 2 23 (18) 

and the sign in the last of equations (|17j) is chosen corresponding to the branch 
of the 2 generation model that one is expanding about i.e A 2 =0 for xi — 
This is because we solved for the coefficients r\,r 2 in terms of square roots of ratios 
of TrA 2 (f = u, d, I). So we should retain an option for the sign of the square root 
: which must be exercised as noted, for consistency, when considering the equations 
appropriate for expansion around the two branches of the toy model solution discussed 
above. The necessity of our choice can be easily understood by considering the two 
generation example where it can be easily checked that A 2 cannot vanish on the (-) 
branch. 

Defining 



A u , 



d,l 



A, 

M 



x 

Sx 



ZZ d + VZ{D T 
3TrZ d + TrZ t 

X - C T -^-C 

TrZ„ 



(19) 



we have finally 



■% = S x + r 2 T ^iS,- {TrZ i,- TrZ ' ) X) (20) 
l rzj u l r/j d 

Since the up sector fermion masses are about 10 2 times larger than the down or 
charged lepton sector it is convenient to rescale the values of the up sector masses 
by 10 2 : TrZ u = 10 2 TrZu and f 2 = 10 2 f2. The definitions are chosen so that when 
S 2 = then S 3 = reduces to Si = 0. Then, in view of the encouraging results of 
our 2 generation model for the 23 dominated fermion hierarchy, it is natural to look 
for solutions to the fitting problem not ab initio, i.e not for a solution of hopelessly 
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non-linear equations but by an expansion (in powers of e ~ \/0~>3 ~ #12 ~ .2) around 
the 2 generation results. Since the 2 generation case gave S 2 ~ 0(e 3 ) it follows that : 

/ = ^r 2 TrZ d S 2 (21) 

The O(10~ 2 ) suppression provided by the ratio d%/v in the above equation implies 
that (/) ~ .01e 3+<5 when f 2 ~ e s . This will then ensure the enhancement of Type I 
neutrino masses that is the rationale for this fitting scenario. 

We chose the SO (10) basis where R = 1 and expand all the orthogonal matrix 
angles 9 u> d,i,T>, (pu,d,i,v, Xu,d,i,v as well as the free parameter f 2 in powers of e. The 
expansion is started by assuming either 

xT = X<°> = Xf ] = X§> (22) 
corresponding to the Xi = +Xi solution, or 

XT = Xl 0) = -x! 0) = ~XS ] (23) 

corresponding to the xi = ~Xi solution. Furthermore, since the off-diagonality 
of the sum rule eqn. (|8j) is driven by that of the CKM matrix which has C ~ e 3 we 
make the weak and well justified assumption (further justified post hoc by actually 
finding solutions of this type) that the angles 4> u ,d,i,v are 0(e) or smaller and thus 
their expansions begin at order e. 

For order e 0,1 ' 2 the expansion simply reproduces the expansion of the toy model 
already described above in Section 3 together with some additional constraints on the 
new angles. For example for the (+) type solutions, at order e°, one gets 9® = 9% 9% = 
9f — 9®. At every order of perturbation we always solve for parameters in which the 
equations at that order are linear and of as high order as possible. This obviates the 
difficulty of multiple solutions of nonlinear equations and their extreme sensitivity 
to the grossly mismatched coefficients implied by the fermion mass hierarchy which 
stretches over 5 order of magnitude. At every order beyond e 7 (where the smallest 
of the mass parameters i.e ui ~ e 7 enters the S3 = equation) we evaluate the x 2 
parameters for the charged fermion fit by setting all undetermined parameters to zero. 
The expansion is pursued till as high as 20th order (!) and the solutions with smallest 
X 2 retained. Then using the truncated values at the chosen order and the values of 
the determined expansion coefficients we reconstruct the fermion mass data via the 
reconstructed Yukawa couplings and mass relation coefficients. We defer a report of 
the details of these solutions to but two points are crucial : the value of x°d an d 
the constraint = I3 — l 2 + d 2 are the same as in the toy model. Thus in line with 
our discussion of that model we shall fit the central values of the charged fermion 
masses and mixing angles except for the value of d 3 which we shall take to be exactly 
d 3 = m r — m M + m s . 
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Till order e 2 the equations are S\ — §2 — 0, thereafter since fi ~ e 3+<5 if ^2 ~ e 5 
it follows that the complete equation S3 = reduces to §x = till one reaches order 
£ 6+<5 Thereafter the full equation S 3 = must be used with the expansion of r-i 
beginning at order e 3+s . Since the up quark mass enters the equations only at order 
e 7 there is no question of finding accurate solutions before that order. Furthermore 
many of the angle expansion coefficients are left undetermined at any given order 
after all the equations available at that order have been solved. These undetermined 
coefficients are simply set to zero before evaluating the results numerically. Note that 
this cannot affect the Orthogonality of the matrices in which the relevant angles occur 
and that this luxury is a consequence of our parametrization of the sum rule. Such a 
procedure may look mathematically questionable particularly since the convergence 
appears to be very weak and typically the values of the masses and mixing angles 
oscillate in an apparently haphazard manner. The reason for this may actually be 
the absence of complexity or the poorness of e as an analytic expansion parameter. 
Nevertheless the proof of the pudding lies in the eating of it ! Our aim -as that of any 
competing numerical fitting procedure - is to find a "close" fit i.e to obtain values of 
the coefficient matrices /, g, h and coefficient quantities r 1; r 2 , r 5 , r 6 , r 7 such that when 
we use the generic form of the SO(10) GUT mass formulae we generate charged and 
neutral fermion masses compatible with the low energy fermion data extrapolated 
over an MSSM desert to the GUT scale. As in the fitting done using the downhill 
simplex method ^3 EH] a simple criterion to pick out an acceptable solution is to 
define \ 2 functions and ask that they be appropriately small : 



2 ■sr^ /T^i TTli \2 

2 v sinflq-sinfl, 2 

XcKM ~ V 5smd a ] {M) 

2 _ 2,2 

Xtot — Xm + XCKM 

For the central quark masses and angles we used the Das-Parida central values 
the fermion data (for tan f3(M s ) = 55 and at M x = 2 x 10 m GeV ) (except for (i 3 as 
explained above) : 



m u = 


.000724 


; rh c = 


.2105 ; 


m t = 


95.148 


rhd = 


.001497 


; rh s = 


.0298 ; 


m b = 


m T + m s — = 1.5835 


m e = 


.000356 


; = 


.0753 


rh T = 


1.629 


sin 9 C = 


.2248 ; 


sin Xc = 


.03278 


; sin (, 


f> c = .00216 



For the 1 — a error estimates we used (all masses in GeV) : 
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= 40.0 ; 


5m c = 


0.018 ; 


5m u = .00013 




= 0.34 ; 


d~m s = 


.0042 ; 


5m d = .0003 


Sm T 


= .038 ; 


5m^ = 


0.00013 


; 5m e = 0.0000015 


5 sin 6 C 


= .0016 


; 5 sin q 


i c = .0005 


; 5 sin Xc = -0013 



For each parameter at an average of (say) half a 1 — a deviation away from 
the central value one might expect a contribution to \ 2 of , roughly .25 so that 
Xm ~ 2.25, Xckm ~ -75- Determining fits much more accurate than these values 
is academic as regards the actual values. It is informative only inasmuch as the 
ability to find finely matched fits to an a priori un-sacrosanct set of central values 
for only approximately known parameters implies only that any such set of data is 
likely to be achievable via a fit of the type indicated. In view of the fact that we 
pushed mj up from its central value by about Srrib/2 to begin with one might even 
expect larger deviations from the "correct" values (the quotes reflect the fact that 
the internal correlations of central values and errors in the extrapolated data at Mx 
are at present very poorly understood). Moreover since the correct scenario for 3 
generations must also explain CP violation, it may be that the slow and oscillatory 
convergence we observe is due to an obstruction in the approach to the true solution 
due to a projection onto a real subspace of the actual parameter space. Surprisingly 
our weakly convergent perturbative method frequently gives solutions with Xtot as 
small as .2 while solutions with x 2 ~ 2 are hardly difficult to find. In the next section 
we discuss examples of the solutions found and of how much closer they have brought 
us to realizing our scenario of a fully successful, NMSGUT compatible, fermion data 
fit. 



3.3 Neutrino Masses and Mixing 

As shown and discussed in |16| |2] the the Type I seesaw mass formula may be 
written 

Ml ^ ^(h + r' 5 g-3f) T r 1 (h + r 5 'g-3f) 



v 

2'a 



= (1.70 x 10-W)iT nRFj sin (3 
= L T VD U V T L (27) 

Where V is the Lepton mixing (PMNS) matrix[2I] in the basis with diagonal leptonic 
charged current and D v the light neutrino masses extrapolated to M X - 

The PMNS matrix V and the ratio of solar to atmospheric mass squared splittings 
can then be identified as 

V = DW ; ^ = &Z^-\ (28) 
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where Af diagonalizes h : n = AfhdiagN' 1 '. Since it is very reasonable to assume 
that the coefficient function Fj will generically have the same magnitudes in the 
NMSGUT as it did in MSGUT i.e < 1, it follows that we can can verify whether 
the enhancement of n makes the NMSGUT generically viable simply by comparing it 
with the typical values n max ~ .5 which failed (for generic Fj values) in the MSGUT 
case by about two orders of magnitude and for special ones by about one order of 
magnitude. 

4 Examples of 3 generation fits 

Our purpose in this contribution is no more than to outline the likely features of 3 
generation fits that follow the scenario suggested [TJ |2] by us as a route to overcoming 
the generic debility of neutrino masses fOJ E] i n the Babu-Mohapatra JU] program. 
Therefore we here only give as examples a few quasi-realistic fits, of the type sought, 
and a few others for comparison, obtained using the method outlined in the previous 
section. As described there we simply took the central values of and, after imple- 
menting the constraint of b = t + s — /i unification ^ required by the consistency 
of the fitting equations, aimed for the central values of the other parameters. There- 
after we fit to the neutrino data as closely as possible. This is what we report here. 
No attempt was made by us as yet to search the parameter space for "optimal" fits. 
A detailed numerical survey will be reported in [21]. Of course improvement of the 
procedure regarding the connection to the low energy data is also possible, but, in 
our opinion this will be called for only after the gross features of the viable fermion 
data fit are established and CP violation is included. In our view this fit should not 
be and is not so unstable and delicately poised as to change its qualitative features 
due to shifts (due to better RG flow control) of a few per cent in the central values 
that one is trying to fit. 

Our procedure was first to fit the ratio of neutrino mass squared splittings to 
the current [22] central value of the solar to atmospheric mass squared splitting ratio 
^ m Li/ ^ m atm ~ -032 by adjusting and thus fixing the value of r' 5 . This gave a number 
of candidate r' 5 values for a given charged fermion data fit. Each value r' 5 corresponds 
to a certain neutrino mass hierarchy and determines a certain value for the eigenvalues 
of n : the crucial matrix whose maximum eigenvalue and coefficient in the Type I 
seesaw mass were generically found^HJ H] to be too small in the MSGUT. Using the 
values of the angles of the orthogonal matrix T> found as part of the solution procedure 
it is easy to generalize the formulae of Section 3 and obtain the three PMNS mixing 
angles. 
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4.1 (+) Type fit 



In the case of the (+) type solutions we find that the solutions obtained assuming 
f 2 ~ g ' 1 ' 2 all yield final f 2 values that differ between the solutions with lowest x 2 by no 
more than about 15%. The x 2 values are also very similar ( x\ot ~ -14 or about .01 per 
variable ! : an excellent fit by normal standards). Since the effect of f 2 on the other 
quantities in the fit is extremely weak, the resultant values of the Yukawa couplings 
and fermion masses and angles are almost identical. This strengthens our notion of 
a unique stable realistic solution/map, of the experimentally found hierarchical type, 
tying together the Yukawa coupling matrices and coefficients to the fermion data. 
Nevertheless a detailed survey is required to determine the actual ranges of mixing 
angles and neutrino mass-squared splitting obtainable and even perhaps alternative 
branches of fits that may actually arise. The values given below have been truncated 
to 4 decimal places but the actual calculation was done retaining 10 digit accuracy. 
In the representative 1+ case where f 2 is assumed to be ~ e to begin with we find 
(fi = (4r 4 )- 1 ) 



15.27 



f 2 = 0.255 



r 6 = 0.0187 ; r 7 = 0.023305 



/ 



/ -0.0002988 0.0000394 
0.0000394 0.009291 
-0.0222691 -0.134823 



V 

/ 0.0000132 
0.0000277 

V -0.0001 

t ° 

-0.00164 

V 0.02478 



0.0000277 
-0.0000146 
0.000057 

0.001644 
0. 

0.119710 



-0.0222691 
-0.134823 
0.481333 

-0.0001 
0.000057 
0.000040 

-0.02478 
-0.119710 




(29) 



(30) 



The eigenvalues of h, f,g are 



h : .518 ; .0277 ; 1.92 x 10" 5 

/ : 1.33 x 10" 4 ; 1.11 x 10~ 4 ; 0.171 x 10~ 4 

g : ±0.122 ; (31) 

Note how the premises of our scenario |2j are indeed respected. 

Then the reconstructed values of the charged fermion masses are(in GeV) 



M v = {95.148,0.2105,0.00077} 
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M D = {1.5835,0.02981,0.0015} 
Mi = {1.629,0.0753,0.000356} 



(32) 



note that we always fit the charged lepton masses exactly. The reconstructed 
CKM matrix is 



C 



0.974449 
0.224548 
0.00526452 0.0324303 



-0.224599 0.002155 
0.973911 -0.032784 
0.99946 



Which yields the CKM angle magnitudes 



(33) 



012 = 0.2265 ; 8 13 = 0.002155 
Then one finds for this solution 



e 23 = 0.03279 



(34) 



Xr 



0.126 



XCKM 



.016 



Xtot 



0.142 



(35) 



To determine the neutrino mixing data the free parameter r' 5 must be fixed. If it 
is known one can calculate n and then since T> is known from the solution the mixing 
angles are easily determined. We fix r' 5 by enquiring for which values of r' 5 the six 
possible ratios 



R 



ijik 



nf -h 2 k 



can attain the current 122 central value R 



m 



sol 



jm 



2 

atm 



(36) 



.32. It is then evident 



that then the indices j,i,k correspond to the e, /i, r neutrinos respectively and so we 
labelled such a solution as an ijik type solution. The software we used allots the 
index 1 to the largest eigenvalue, 2 the next largest and 3 the smallest, this accounts 
for the solution labels encountered below. 

For the I+case in question one finds the following set of solutions : 



J + 2321 : 


(a) 


i 


= -0.33285 


(b) r' 5 


= 0.48265 


7 + 3231 : 


(a) 


i 

r 5 


= -0.33140 


(b) r' 5 


= 0.481314 


J + 2123 : 


(a) 


i 

r 5 


= -0.14427 


(b) 4 


= -0.13932 




(c) 


i 

r 5 


= 0.301785 ; 


id) r' 5 = 


= 0.30650 


7+1213 : 


(a) 




= -0.14435 


(b) r' b - 


= -0.139237 




(c) 


r' 5 


= 0.306578 ; 


id) r> 5 


= 0.30171 



It is clear that the cases 1+2321 and 1+3231 and cases 1+2123 and 1+1213 are 
almost identical. Note also how close I+2123(a) and I+2123(b) and I+2123(c) and 
12123(d) are. In fact these pairs give almost identical leptonic mixing parameters 
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and n eigenvalues. Case 1+2321 corresponds to the hierarchy m Vr » m v » m Ve 
while case 1+3231 corresponds to m„ r >> m Ue » m v . Similarly case 1+2123 
corresponds to the hierarchy m Ue » m Ufi » m VT and case 1+1213 to m VfL » 
m Ve » m VT . Cases 1+2321,1+3231 have highly hierarchical n eigenvalues while in 
cases 1+2123,1+1213 which have an inverted hierarchy there are two large almost 
degenerate eigenvalues and one very small one. Because of the close similarities we 
give only the values for cases I+2321a,I+2321b and I+2123a,I+2123c in Table 2. 
Only in case 1+2312 is the 23 sector mixing angle large and the 13 mixing angle 
reasonably small. However the value of the 12 sector mixing is very low. Clearly then 
two large and one small PMNS mixing angles are hard to achieve. On the other hand 
the large ( about 200 times larger than the Type I fits in the 10 — 126 scenario) value 
of the largest n eigenvalue together with the satisfactory value of the mass squared 
splitting ratio means that the problem with too small neutrino masses is unlikely to 
appear even for generic values of the GUT scale breaking, leave alone regions where 
the coefficient function Fi is itself large. 

Even though they are "captured" at different orders in perturbation theory the 
solutions obtained when f 2 is assumed to be ~ 1 or ~ e 2 to begin with are very close 
to the one displayed above which was obtained by assuming that f 2 ~ e. In fact the 
final values of f 2 in these two cases are .00185, .00225. Some of the details can be 
compared in Table I. 

Focussing on the case 1+23 12 (a) we find that using the value r' 5 = —.333 and the 
values of the T> angles found for the 1+ solution : 

<fo = 0.0525 ; 9 V = 0.43481 ; xv = 0.02568 (38) 
then gives, on using eqn. (j28|) the following values for the eigenvalues of h 

hi = 113.8 ; n 2 = 20.045 ; n 3 = 1.97 x 10~ 5 (39) 

Which clearly shows the required 100-300 fold enhancement relative to the BM 
fitting scenario in the overall scale of the Type I seesaw masses via an enhancement 
of h that motivates our scenario. However the PMNS angles are found to be : 

sin 2 #f 2 M7VS = .073 ; sin 2 fl 2 p 3 M7VS = .77 ; 13 = .176 (40) 

Although one might hope to improve the value of the 23 and 13 angles towards 
more realistic values it is clear that raising the 12 sector mixing from such a small 
value will be difficult if the fit is indeed stable as we propose it is. A final conclusion 
must await a detailed survey. However recall that even in the 10 — 126 scenario 
no Type I solution was found fTJ till CP violation was included in the charged 
fermion sector. Results of the lepton sector parameters for some for some of the other 
distinct cases are given for comparison in Table 2. 
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1+ 


11+ 


III+ 


I- 


II- 


Hi- 


N 


15 


15 


9 


18 


20 


ll 


5 


1 





2 


1 





-1 


f 2 


.255 


.225 


.185 


1.21 


1.065 


- 5.07 


r e 


.01867 


.01867 


.018667 


.01886 


.01886 


.018913 


r 7 


.0233 


.0233 


.0233 


-.0234 


-.0234 


-.0233 


k 


.518 


.518 


.518 


.519 


.519 


.519 


k 


.0277 


.0277 


.0277 


.02689 


.02687 


.02651 


10 5 /l3 


1.92 


1.99 


2.05 


16 


16.55 


24.5 


9 


±.122,0 


±.122,0 


±.122,0 


±.120,0 


±.120,0 


±.120,0 


10 4 /i 


1.33 


1.17 


.961 


5.39 


4.705 


16.67 


io 4 / 2 


1.11 


0.98 


0.80 


4.70 


4.10 


14.85 


io 4 / 3 


0.171 


0.150 


0.12 


0.76 


0.66 


1.60 




0.126 


0.148 


0.128 


0.06 


1.56 


0.255 


XCKM 


0.016 


0.008 


0.0132 


1.641 


0.075 


0.051 



Table 1: Charged fermion fitting parameters and Yukawa eigenvalues for 6 sample 
fits. N is the order of perturbation theory at which the fit was "captured" while the 
exponent 5 specifies the initial assumption f 2 ~ e s . 
To the accuracy displayed the value of f x = 15.27 for every solution. 

4.2 (-) Type fits 

For the branch defined by imposing Xd = Xd = ~ X? = ~Xv on the perturbative 
iteration, the fits we obtain are, at \ 2 > 1-6, somewhat poorer than the ones found 
for the (+) case. Moreover even such fits first occur at quite high order in the 
perturbation expansion and the differences between the final f 2 values are somewhat 
greater. The acceptable cases with final values f 2 ~ .01 are shown in Table 1. as cases 
I- , 1 1- along with another fit III- where an initial value f 2 ~ e _1 ~ leads to a distinct 
solution with f 2 ~ .05. This type of solution has n eigenvalues {4.3, .75, .012} and 
shows clearly that / < .Ole 3 is necessary to raise the value of n enough to overcome 
the difficulty uncovered for the MSGUTjTBllI]. The rest of the procedure is identical 
to the (+) case and the relevant results are collected in Tables 1 and 2. 

For the readers convenience we give the values of the Yukawa couplings, coefficients 
ti for one of the (-) cases. Namely I-2312a. 



15.27 



f 2 = 1.21 



r 6 = .01886 



r 7 



-.02340 



/ -0.000149706 -0.00054677 -0.014659 
-0.00054677 0.00927299 -0.134398 
V -0.014659 -0.134398 0.48307 
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I+2321a 


I+2321b 


I+2123a 


I+2123c 


Hierarchy 


t > /i > e 


t > /i > e 


e > /i > t 


e > fi > t 




-.333 


.483 


-.144 


.302 


hi 


113.83 


104.26 


54.38 


51.95 


h 2 


20.045 


18.36 


52.53 


51.14 




1.97 x 10~ 5 


1.89 x 10~ 5 


1.78 x 10~ 5 


1.774 x 10- 5 


sin 2 2# 23 


.768 


.553 


.0726 


.077 


6*13 


.176 


.117 


1.347 


1.35 


sin 2 2^i2 


.073 


.130 


.143 


.753 




I-2312a 


I-2312b 


I-1223a 


I-1223c 


Hierarchy 


t > /i > e 


t > /i > e 


e > /i > t 


e > /i > t 




-.258 


.698 


-.073 


.569 


hi 


32.81 


20.96 


14.99 


11.12 




5.78 


3.69 


14.75 


10.95 


h 3 


4.65 x 10" 4 


3.36 x 10" 4 


4.53 x 10~ 4 


3.76 x 10" 4 


sin 2 2# 23 


0.886 


0.447 


0.192 


0.175 


#13 


0.216 


0.159 


1.261 


1.282 


sin 2 2# 12 


0.209 


0.206 


0.494 


0.772 



Table 2: Values of neutrino mass and mixing parameters determined assuming central 
value for the ratio of solar to atmospheric mass splitting. 

I 0.0000521467 0.000139 -0.00044 \ 

0.000139 -0.000039 0.000182 (41) 

V -0.00044 0.000182 0.000133 / 

/ 0.00116 -0.01753 \ 
-0.00116 -0.1191 

V 0.01753 0.1191 / 
-0.04057 ; 6 V = -0.134924 ; X v = -0.57541 

Mu = {95.148,0.2105,0.00077} 
M D = {1.5835,0.0298017,0.0014962} 
Mi = {1.629,0.0753,0.000356} (42) 

5 Discussion, Conclusions and Outlook 

In this paper we have reported progress in accomplishing our program EEl H 12] 
for a completely realistic fit of all known charged fermion and neutrino mass data 



/ = 
9 = 

4>v = 
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using the mass relations and RG evolution common to any SO(10) Susy GUT with a 
10 — 120 — 126 FM Higgs system fermion. Specifically, we have shown that in the 
quasi realistic 3 generation but CP preserving (real) case we are able to obtain accurate 
charged fermion fits, neutrino mass parameters and a PMNS mixing pattern that can 
be large at least in the 23 sector and small in the 13 sector. The remaining deficiencies, 
namely the absence of simultaneous large mixing in the 12 sector and the fit of the 
MSSM CKM CP phase in the first quadrant can presumably be remedied in the 
complex 3 generation case, in close analogy with the BM program where a succesful 
Type I fit could be found [TT) ITo] only when CP violation was introduced. For the case 
where the couplings are all real and CP violation arises from the phases present in the 
coefficients the CP violation can be traced to the GUT scale symmetry breaking 
and fine tuning^ HH Hj that defines the light Higgs pair of the MSSM. In this case 
the EW doublet ([1, 2, ±1]) mass matrix formulae for the NMSGUT that we have 
already furnished will permit the evaluation of the feasibility of the generic fit for 
realizability in the parameter space of the NMSGUT. The report on the investigation 
of these questions will be given in [26J. 

We carried a somewhat novel expansion in a single parameter e, in terms of 
which the very strong (~ 10 5 ) hierarchy visible in the MSSM fermion data can be 
compactly formulated, to high orders of the perturbation expansion . The explicit 
and systematic incorporation of the hierarchical structure into our fitting procedure 
makes it well suited to verification of the viability of the hierarchical fitting sce- 
nario we proposed^ Ej. On the other hand the convergence in the parameter e 
in the three generation case(in sharp contrast to the 23 generation toy model[2|) is 
only oscillatory |24j. The 2 generation Toy model is necessarily CP conserving in the 
charged fermion sector whereas the true solution in the 3 generation case is not. Thus 
the failure to converge may reflect an obstruction to fitting the data to the absolute 
minimum when the channels to reach it on the complex plane have been blocked 
by our reality assumption. Our method is complementary to alternative numerical 
procedures such as the "downhill simplex" method[lfl| ITH] which it can verify as well 
as provide initial search points for. Inasmuch as it is only the actual hierarchical 
mass data that must be fitted the true solution must necessarily be "close" to the 
23 generation toy model caricature. Thus the "missed solutions anxiety" associated 
with these methods, may find some palliation when the fitting problem is considered 
in the light of the insight provided by these two complementary methods jointly. It 
is only after the successful completion of this program that we may, finally, be in 
possession of a well defined Susy GUT compatible with all low energy data as well 
as information on the embedding of the MSSM in the MSGUT : which emerges as 
the most valuable corollary product of the fitting procedure JH]- It is only then that 
we will be able to enter the third phase of the GUT program in which the exotic 
process (AB ^ etc ) predictions will finally begin to be formulated in a sufficiently 
tight manner as to make the comparison with data a falsifiability test rather than a 
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hopeful check on a lottery bet. 
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7 Note Added 

On 18 July 2006, as this paper was being readied for submission to the archive, we 
received |2Zj in which a successful fit to the complete fermion data including CP 
violation, for tan/3 ~ 10, with spontaneous CP violation, and with an additional Z 2 
symmetry imposed to further reduce the number of degrees of freedom, is reported. 
The results are consonant the scenario proposed in [IJ E] and the results of this paper. 
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